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ABSTRACT 


Context. Transitional disks with large dust cavities are important laboratories to study planet formation and disk evolution. Cold gas 
may still be present inside these cavities, but the quantification of this gas is challenging. The gas content is important to constrain the 
origin of the dust cavity. 

Aims. We use Atacama Large Millimeter/submillimeter An'ay (ALMA) observations of '^CO 6-5 and 690 GHz (Band 9) continuum 
of five well-studied transitional disks. In addition, we analyze previously published Band 7 observations of a disk in '^CO 3-2 line 
and 345 GHz continuum. The observations are used to set constraints on the gas and dust surface density profiles, in particular the 
drop dgas of the gas density inside the dust cavity. 

Methods. The physical-chemical modeling code DALI is used to analyze the gas and dust images simultaneously. We model SR2L 
HD135344B, LkCal5, SR24S and RXJ1615-3255 (Band 9) and J1604-2130 (Band 7). The SED and continuum visibility curve 
constrain the dust surface density. Subsequently, the same model is used to calculate the '^CO emission, which is compared with the 
observations through spectra and intensity cuts. The amount of gas inside the cavity is quantified by varying the (5gas parameter. 
Results. Model fits to the dust and gas indicate that gas is still present inside the dust cavity for all disks but at a reduced level. The 
gas surface density drops inside the cavity by at least a factor 10, whereas the dust density drops by at least a factor 1000. Disk masses 
are comparable with previous estimates from the literature, cavity radii are found to be smaller than in the 345 GHz SubMillimeter 
Array (SMA) data. 

Conclusions. The derived gas surface density profiles suggest clearing of the cavity by one or more companions in all cases, trapping 
the millimeter-sized dust at the edge of the cavity. 
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1. Introduction 

Planets ar e formed in disks of dus t and gas that surround young 
stars (e.g. I Williams & Ciezall20li1) . Of particular interest are the 
so-called transitional disks with cleared out cavities in the inner 
part of the dust disk. These disks are expected to be in the middle 
of active evolution and possibly planet formation. Transitional 
disks were originally identified through a dip in the mid in¬ 
frared part of the Spectral Energy Distribution (SED), indicat¬ 
ing a lack of small micron-sized dust grains close to the star. 
They wer e modeled as axisymmetric dust disks with inner cav¬ 
ities (e.g. IStrom et ^Il989l iBrown et"^ 120071: lEsnaillat et al.l 
l2014l) subsequently confirmed by millimeter interferometry 

(e.g.lPu trey et al .l2008l:lBrown et al.l2009l: [Andrews et al.l201 It 

It _ _ n _ _ _ I i ^-tK t.. __ i._i. . . _ .^i. _ . . _ 


llsella et al.l2012h . In recent years it has become clear that the gas 


itv (ISalvk et a 

l.l2007t Pontoppida 

n et a1.l2008HSa1vk et a1.l2009l 

iBrittain et al.l 

2009i IBrown et al.l 

2012h. However, these data do 


tional line obser vations reveal directly the presence of gas inside 
the dust cavity ([van der Marel et al. 20131: Bruderer et Sl l20l4 


iPerez et ^l2()14 : ICasassus et al.l 201 31: Zhang et al .1120141) . Also 

the size distribution of the dust is possibly not uniform through¬ 
out the disk: micron- and millimeter-siz e d grains may have 
differ ent distributions (iGarufi et al.l 1201 3t Ivan der Marel et akl 
120131) and mill imeter- and centim eter-sized grains can show ra¬ 
dial variations (iPerez et al.ll2012l) . 

Quantifying the density structure of the gas inside the dust 
hole is essential to distinguish between the possible clearing 
mec hanisms of the dust in the inner par t of disk: photoevapora¬ 
tion dClarke et al.l 1200 ll) . grain growth (iDullemond & DominikI 


20051), clearing by (sub )stellar or planetary companions (e.g. 


Lin & Papaloizoull 19791) . Also, instabilities at the edges of dead 


and dust do not neces sarily follow the same ra dial distribution 
in transition disks (e.g. lPont oppidan et al.ll200 8). High accretion 
rates (~ 10 * Mq yr"*, e.g. Calvet et al.l2005l: Naiita et al.ll2007l) 


indicate that the dust cavities are not empty, but gas is still flow¬ 
ing through the gap. Near infrared observations of CO rovibra- 


not trace the bulk of the molecular gas. ALMA CO pure rota- 


Jansky Eellow of the National Radio Astronomy Observatory 


zones ha ve been involved to millimeter-dust rings and asymme¬ 
tries (e.g. lRegalv et al.l2012l) . Planet cand idates have been foun d 
in cavitie s of the transition disks T Cha (iHuelamo et^ 20nE 
LkCa 15 (iKraus & Irelandll2012l) and HD 100546 (lOuanz et al.l 
l2014l) . indicating planet clearing as a possible scenario. 

Since detection of planets in disks is challenging, an indica¬ 
tor of the clearing by a (planetary) companion can be inferred 
from the different structure of the gas compared to the dust in 
several of the observed disks as shown by planet-disk interac- 
tion models in co mbination with dust evolution models (e.g. 
iPinilla et alJl2012l) . Companion clearing changes the gas surface 
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density profile from that of a full disk; the density at and around 
the orbital radius of the planet is decreased while at the outer 
edge of this gas gap the density slightly increases. Due to the in¬ 
duced pressure gradients along the sides of this pressure bump, 
dust particles will drift towards the edge of the gap, get trapped 
and grow efficiently to larger sizes. This r esults in a dust rin g 
of millimeter-sized dust or radial dust trap (iPinilla et al.ll2012]i . 
The decrease in gas density due to companion clearing depends 
on the mass of the companion and other parameters such as 
the viscosity (IZhu et al.ll201 lUPodson-Robinson & S alvkll201 ll 

iPinilla et al.ll2012l : lMulders et al.ll2013HFung et al.ll20l3r Other 

effects that play a role in dynamic clearing a re for example 
migra tion and planet accretion luminosity (e.g. iKlev & NelsonI 
l2012h . Therefore, the amount of gas inside a cavity can help to 
constrain the mass of the companion in the planet disk clearing 
and dust trapping scenario if the viscosity is known. In compar¬ 
ison, the surface density of the millimeter-sized dust drops sig- 
nihcantly deeper inside the cavity compared to the g as due to the 
trapping of millimeter-sized dust (IPinilla et al.ll2012l) . Therefore, 
the quantihcation of gas and dust inside the cavity indicates if 
the trapping mechanism due to a planet is at work. Dust trapping 
at the edges of dead zones will also result in millimeter-sized 
dust rings, but in that scenario the gas surface density remains 
the same as for a full disk, apart from a min or density increase 
(less than 50%) at the edge (iLvra et alJl2015li . 

The rovibrational lines used in previous work of transitional 
disks can indicate the presence of CO inside the dust cavity 
and the inner radius, but are not suitable for quantihcation of 
the total gas density. The emission is very sensitive to the gas 
temperature, which is high but uncertain in the models in this 
regime. Also, the dust opacity of the inner disk, which provides 
shielding, is uncertain. The pure rotational CO lines provide bet¬ 
ter constraints on the gas surface density inside the cavity. The 
brightest molecular lines in the submillimeter regime are the ro¬ 
tational *^CO lines. Before ALMA, interferometric data gen¬ 
erally did not have the sensitivity and spatial resolution to de¬ 
tect rotational '^CO emission inside the cavity. However, even 
if CO line emission is detected, this emission can not be di¬ 
rectly translated into densities, especially inside the dust cav¬ 
ity: the CO abundance is not constant throughout the disk with 
respect to molecular hydrogen due to photodissociation by the 
UV radiation from the star and free ze out onto dust grains in 
colder regions clos e to the midplane (Ivan Zadelhoff et ani200lt 
lAikawa et al.ll2003) . Furthermore, the gas temperature is decou¬ 
pled from the dust temperature in the surface of the disk, particu¬ 
larly at the cavity wa ll that is directly heated and inside the dust- 
depleted cavity (e.g. Kam p & Dullemondll2004t iJonkheid et alJ 


, _ imp ( 

l200llGorti & HollenhachlboOSl) . 

To determine the amount of gas that is present inside a dust 
cavity, the physical and chemical structure of the gas needs to be 
modeled. For this we make use of the physical-chemic al mod¬ 
eling tool DALI (iBruderer et al J l2()T^ iBrudereil l20 1 3l) . which 
solves the heating-cooling balance of the gas and chemistry si¬ 
multaneously to determine the gas temperature, molecular abun¬ 
dances and molecular excitation for a given density structure. 
Due to the disk Keplerian rotation and gas thermal and turbulent 
motions, the emission along any line of sight that inter¬ 
cepts the disk is a spectral line comprising of an optically thick 
central part and potentially optically thin wings. With spatially 
resolved CO observations it is possible to infer emission from 
different line s of sight, in contrast to previous single dish spec¬ 
tral analysis. IBruderer et al.l (l2014li demonstrated for 6-5 
lines in Oph IRS 48 that the emission in the optically thin wings 
can be used to constrain the gas density. In this case the observa¬ 


tions are found to be consistent with a physical disk model with 
two gas density drops inside the millimeter dust cavity. However, 
Oph IRS 48 has a very low disk mass (constrained by C'^O ob¬ 
servations) so even the ^^CO line is close to optically thin. For 
more massive disks, the bulk of the *^CO emission may remain 
optically thick even if the gas density inside the dust depleted 
cavities drops by one or two orders of magnitude. In that case, 
the '^CO emission originates from deeper in the disk (closer to 
the disk midplane) where the gas is colder then in th e surface, 
causi ng a drop in the optically thick line intensity (iBrudereil 

l2()H . 


To date, a total of 9 transitional disks with large resolved 
millimeter dust cavities has been observed in ALMA Cycle 0 
in dust continuum and '^CO. In this paper, we present observa¬ 
tions of the disks SR21, HDI35344B, RXJI 6I5-3255, LkCal5 
and SR24S through '^CO 6-5 observations (iPerez et al. 2014ll 
and J1604-2130 through '^CO 3-2 observations ( Carpenter et alJ 
120141 : IZhang et~^l2014l) . With DALI we determine the density 
structure of the gas and dust that is consistent with the SED, dust 
continuum interferometry and the '^CO observations. If avail¬ 
able, information on the hot gas and dust from the literature is 
included. The goals of this study are to determine the amount 
of the gas and dust inside the dust hole and the implications for 
potentially embedded planets, assuming that the millimeter dust 
rings are indeed dust traps. The derived physical models in this 
study can be used as a starting point for future ALMA observa¬ 
tions of CO isotopologues, whose emission is optically thin and 
which therefore directly trace the density structure. Our analysis 
of the gas density structure is similar to that used for Oph IRS 48 
th rough *^CO 6-5 and C ^^O 6-5 line observations as presented 
in lBru^derer et alJ (l2014li . The analysis of HD142527, which has 
been observed in ^ ^CO, ^^CO and C^^ O 6-5, 3-2 and 2-1, is left 
for future studies. iPerez et ^ (1201 5l) have modeled the 7=2-1 
CO isotopologues and determined the gas mass remaining inside 
the cavity of that disk. Other transition disk s that have been ob - 


served in ALMA Cycle 0 are H D 100546jWffi^h_etffi 


TW Hya (iRosenfeld et al.ll20l4 . Sz 91 (E anovas et ^ l2nL5h 
and MWC 758 (not yet published), but since their dust holes 
remain unresolved in the current data they are not considered in 
this study. 

The paper is structured as follows. In Section 2.1 we describe 
the details of the ALMA observations. In Section 2.2 we present 
moment maps of the *^CO observations and the derived stellar 
position and orientation. The modeling approach and constraints 
from the literature are presented in Section 3. Section 4 presents 
the modeling results. Section 5 discusses the implications for 
embedded planets in the disk. 


2. Data 

The continuum and *^CO line observations were obtained dur¬ 
ing ALMA Cycle 0 between June 2012 and January 2013, with 
baselines ranging from 20 to 500 m, probing scales from 0.2 to 
4 arcseconds. The sources and their properties are summarized 
in Table 121 

2.1. Observational details 

Five transitional disks were observed in ALMA program 
2011.0.00724. S (PI Perez), two of these disks were presented 
in IPerez et al.l (l2014l) . The disks were observed in Band 9 (690 
GHz, 440 yum), with four adjacent spectral windows of 1875 
MHz, each with 3840 channels of 488 kHz (0.2 km s“*) width. 
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Table 1. Stellar properties 


Target 

SpT 

L, 

(Lo) 

M, 

(Mo) 

R. 

(Rg) 

Ttfs 

(K) 

M 

(MoS-‘) 

d 

(pc) 

Ay 

(mag) 

Ref. 

SR21 

G3 

10 

1.0 

3.2 

5830 

< 1 • 10-" 

120 

6.3 

1,2,3 

HD135344B 

F4 

7.8 

1.6 

2.2 

6590 

6 • 10-8 

140 

0.3 

4,5,6 

LkCal5 

K3 

1.2 

1.0 

1.7 

4730 

2 ■ 10-8 

140 

1.5 

6,7,8 

RXJ1615-3255 

K5 

1.3 

1.1 

2.0 

4400 

4-10-'° 

185 

0.4 

5,9 

SR24S 

K2 

2.5 

0.7 

2.7 

4170 

7 • 10-* 

120 

7.0 

5,10,11 

J1604-2130 

K2 

0.76 

1.0 

1.4 

4550 

1■10-" 

145 

0.66 

12-14 


p iDunkin et all dl^ ), 2) Pontoppid an et all 1200^ . 3) iGrady e t all 120091). 4 1 iPrato et all (I20031). 51 lAndrews et all ll201lh. 6) lEspaillat et al 
( 2010 h . 71 lPietu etal.l(l^ ), Slllnglebv eTnToom 91IWichmann et all ( Il997ll . lOl lNatta et alj J2006h ! 1 li lLuhman & Ri^ jl999h . 121 10" 


1 2008 i. 13l lMathews et all 


)■ Sllpateby et al.iazuu'^). tJ)IWich 

20131 . 141 icarpenter et alJ2014ll 


Table 2. Properties of the ALMA observations 


Target 

Derived position 
(J2000) 

Line 

Beam 
size (”) 

Beam 

PA(°) 

rmsii„e° 

(mJy) 

Pcakjjne 

(mJy) 

rmSco„, 

(mJy) 

Peakcoiit 

(mJy) 

PA 

(°) 

i 

n 

^src 

(km/s) 

SR21 

16:27:10.27 -24:19:13.03 

'"CO 6-5 

0.28x0.19 

-85 

65 

1660 

1 

409 

14 

15 

2.75 

HD135344B 

15:15:48.43 -37:09:16.36 

'"CO 6-5 

0.26x0.19 

49 

114 

2010 

1.8 

318 

63 

20 

7.05 

LkCal5 

04:39:17.80 -h22:21:03.21 

'"CO 6-5 

0.30x0.23 

-33 

120 

1202 

1.8 

118 

60 

55 

6.1 

RXJ1615-3255 

16:15:20.23 -32:55:05.40 

'"CO 6-5 

0.29x0.20 

-74 

102 

960 

1.7 

123 

153 

45 

4.6 

SR24S 

16:26:58.51 -24:45:37.14 

'"CO 6-5 

0.40x0.20 

-75 

96 

2196 

1.9 

278 

20 

45 

4.75 

J1604-2130 

16:04:21.65 -21:30:28.90 

'"CO 3-2 

0.69x0.44 

-78 

100 

1500 

0.20 

38 

80 

10 

4.7 


Notes. Measured in 0.5 km s * bins. 


for a total continuum bandwidth of 7.5 GHz. One spectral win¬ 
dow was centered on the '^CO 6-5 line (691.47308 GHz). The 
bandpass was calibrated using 3C279, Titan was used as flux cal¬ 
ibrator and the phase and amplitude were calibrated using J1427- 
4206 (HD135344B), J1625-2527 (SR21, SR24S and RXJ1615- 
3255) and J0510-I-180 (LkCal5). Each target was observed for 
~25 minutes on source. Given the high signal to noise ratio of 
the continuum, amplitude and phase self-calibration were per¬ 
formed on all data and applied afterwards to the '^CO data as 
well. 

Disk J1604-2130 was observed as part of ALMA pro¬ 
gram 201L0.00526.S (PI Carpenter), in Band 7 (345 GHz, 880 
pm), with one spectral window centered on the '^CO 3-2 line 
(345.79599 GHz) with 488 kHz (0.42 km s^^) ch annel width 
The da ta and observational details are presented in IZhang et alJ 
(l2014l) . Due to the longer wavelength, the beam size of these 
data is considerably larger than that of the Band 9 data. 

The data were imaged using CASA v4.2, adopting the Briggs 
weighting scheme with robust of 0.5. The resulting rms and 
beam sizes of the line data are given in Table |2 During the 
cleaning process of the SR24S data, it became clear that there 
is extended emission that could not be recovered by the shortest 
baselines of these ALMA observations (50 kd, or 4” scales). 

2.2. Continuum and line maps 

Figure [T] shows the ALMA continuum image, zero-moment 
'^CO map and '^CO spectrum of each of our targets. The spec¬ 
trum was extracted from the region of the zero-moment map 
size. The contours indicate the 3cr detection limit in each im¬ 
age. The continuum images reveal ring-like structures for all tar¬ 
gets except RXJ1615-3255, even though its SE P and the 345 
GHz visibility curve show signs of a dust cavity (lAndrews et al.l 
1201 ill . The rings clearly show azimuthal asymmetries, which 
have been confirm ed to be real inten sity asymmetries for SR21 
and HD135344B (iPerez et al.llMl . For SR24S and LkCal5 
the asymmetries are likely caused by the geometry of the disk 


{i ~ 45°) and the dust continuum becoming optically thick at 
this wavelength (0.45 mm). The S/N on the continuum varies 
between ~65 and 400, much higher than previous SubMillimeter 
Array (SMA) observations. 

The zero-moment '^CO maps reveal the presence of gas in¬ 
side the dust hole for all targets, demonstrating different distribu¬ 
tions of gas and dust. This does not mean that there can be no de¬ 
crease in the gas density compared to a full disk density profile, 
because the '^CO emission is likely optically thick. The maps 
of HD135344B and J1604-2130 show a decrement of emission 
around the stellar position, directly hinting at a strong decrease 
of gas density. The size of this gas ‘hole’ is much smaller than 
the dust hole and comparable to the beam size. For LkCal5 there 
appears to be a slight offset between the continuum and inte¬ 
grated CO emission peaks; the maps are currently aligned to the 
same central positi on. This may be relate d to the proposed ec¬ 
centric inner disk (iThalmann et al.l [20141) . casting a more pro¬ 
nounced shadow on one side of the disk. However, considering 
the limited significance of this shift, modeling this asymmetry is 
beyond the scope of this paper. 

All *^CO spectra are consistent with a rotating disk at 
an inclination >15°, revealing a double-peaked velocity pro- 
file, except J1604-2130, which is known to be almost face-on 
(iMathews et al.ll2()T3 IZhang et al.ll2014l) . The spectra of SR21 
and SR24S are likely affected by foreground absorption, com- 
monly observed in the Ophiuchus star forming region (e.g. 
Ivan Kemnen et al.l l2009h . although no single dish spectra of 
these targets are available to confirm this. The SR21 spectrum 
is somewhat asymmetric: the red side is narrower than the blue 
side. Furthermore, the data of SR24S are affected by extended 
emission outside the field of view of the ALMA observations. 
This may be related to the suggested circumbina ry disk around 
SR24N, which is located 5” north, and SR24S (iMavama et al.l 
1201 Ol) . but this can not be confirmed with the available data. The 
apparent double peak in the spectrum is caused by foreground 
absorption: the bulk of this emission is not originating from the 
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Fig. 1. ALMA observations of the continuum and '^CO line. The first 5 disks show the 690 GHz (440 continuum and the '^CO 
J—6-5 line, the 6th is the 345 GHz (880 //m) continuum and J—3-2 line. Left: Continuum image. The stellar position is 
indicated by a white star, the white bar in the upper right corner indicates the 30 AU scale and the yellow contour gives the 3cr 
detection limit. The colorbar units are given in Jy beam *. ; Center: zero-moment '^CO map. The colorbar units are given in Jy 
beam^' km s Right: '^CO spectrum integrated over the entire disk. The dashed line indicates the zero flux level and the grey 
areas indicate the parts of the spectrum affected by foreground absorption (seen in SR21 and SR24S). The beam is indicated in each 
map by a white ellipse in the lower left corner. 
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disk. Therefore, we do not aim to model the gas surface density 
of this target and only model the dust continuum structure. 

The '^CO channel maps (binned to 1 km s ') and first mo¬ 
ment maps are presented in Figure lA.ll We derive the stellar 
position, inclination, position angle and source velocity from 
the channel maps and first moment map and overlay the cor¬ 
responding Keplerian velocity pattern in white contours. The 
derived parameters are consistent within 5° with values from 
the literature for resolved gas and dust , wher e available (e.g. 
iPontoppidan et alJ I2008L [Andrews et al.l I 20 TT!) . For the stellar 
mass, we use values from the literature (see Table[TJ. The derived 
parameters are given in Table |2l The peak S /N of the channel 
maps is at least 10. 

The channel maps confirm the presence of gas inside the dust 
cavities. The velocity range that falls within the derived dust cav¬ 
ity radius is given in Table lA.ll The velocity is calculated 

using Vobs = sin / with G the gravitational constant. The 

\ f cav 

channel maps also reveal that the detected CO emission is more 
extended in the outer disk than the dust continuum for all sources 
except SR21, although the signal to noise and image quality are 
low in the outer parts of the disks. The emission is unlikely to be 
affected by spatial filtering, as shown in the visibilities (Figure 
H; the emission is well covered down to ~40k/l, corresponding 
to 5” (>600 AU) diameter scales. Whether or not the difference 
in extent is significant needs to be quantified with a physical 
model, because the dust emission may simply be too weak in 
the outer disk to be detected while '^CO remains detectable be¬ 
cause it is optically thick. The red-shifted side of both SR21 and 
HD135344B is somewhat weaker than the blue-shifted side, as 
also seen in the spectrum (Fig.[Tli. The location of the decrease is 
cospatial with the peak of the asymmetry in the continuum map. 
This is possibly related to a higher continuum optical depth or 
lower temperature at the location of the asy mmetry, as was also 
seen for H 2 CO emission of Oph IRS 48 (Ivan der Marel et aH 
12014 . 

3. Method 

3.1. Modeling 

Our aim is to constrain a physical disk model that describes both 
the dust and gas structure of each disk. The main goal is to de¬ 
termine the drop in gas surface density bg^s inside the dust hole, 
within the constraints of a model that fits the gas and dust in the 
outer disk (Fig. |2]i. Therefore the other disk parameters need to 
be constrained first, as they have a significant influence on the 
strength of the CO emission inside the cavity: the disk mass re¬ 
lates to the optical depth of the ^^CO lines; the vertical structure 
affects the temperature by shadowing and the amount of exposed 
disk surface atmosphere by stellar light; the amount of dust in¬ 
side the cavity, including the inner disk, shields the CO from 
photodissociation but also provides cooling; and the location of 
the cavity wall determines the temperature of the wall, where the 
bulk of the '^CO emission is expected to originate (see 3.3.3 for 
details). To this end, information is combined from the SED, the 
ALMA continuum and integrated CO visibilities, and the CO 
disk integrated spectrum and zero-moment map. 

Although the model is not expected to be unique, it sets a 
structure for the disk that can be tested with future ALMA ob¬ 
servations of CO isotopologues. The main goal of this study is 
to compare the maximum possible decrease of the gas surface 
density with the minimum decrease of the dust density inside 
the cavity. If the drop in gas density is less deep than for the 


millimeter-sized dust, this is direct evidence for the trapping sce¬ 
nario. 

The axisym metric physical- chemical model DALI is used 
for the analysis (lBrudereill2013L details next section) , thus any 
azimuthal asymmetries are ignored. The most prominent asym¬ 
metries are those seen in t he millimeter-size d dust continuum 
in SR21 and HD135344B (iPerez et al.ll20l4 . but the influence 
of the millimeter dust on the gas temperature and UV shielding 
is negligible so axisymmetry can be safely assumed. The small 
dust grains are the main heating agents of the gas. In the visibili¬ 
ties we therefore use the amplitude ( y[{Reaf + because 

the asymmetries and thus the variations in the imaginary part are 
small. The visibilities constrain the radial structure of the emis¬ 
sion. Due to limited u, v-coverage in these Cycle 0 data the vis¬ 
ibilities give better constraints on the structure than the images 
alone, especially on the cavity radius through the location of the 
first null. 


3.2. DALI 

The determination of the gas disk from the '^CO emission re¬ 
quires a thermo-chemical disk model, in which the heating¬ 
cooling balance of the gas and chemistry are solved simulta¬ 
neously to determine the gas temperature and molecular abun¬ 
dances at each position in the disk. Moreover, even though the 
densities in disks are high, the excitation of the rotational lev¬ 
els may not be in local thermodynamic equilibrium (LTE), and 
the gas and dust temperature are decoup led, especially inside 
and at the cavit y edge. The DALI model (iBruderer et al.luOia 
lBrudereJl2013h uses a combination of a stellar radiation field 
with a disk density distribution as input. DALI solves for the dust 
temperatures through continuum radiative transfer from UV to 
millimeter wavelengths and calculates the chemical abundances, 
the molecular excitation, and the thermal balance of the gas. It 
was developed for the analysis of gas emission structures in pro¬ 
toplanetary disksJiKludi^ transitional disks with varying gas to 
dust ratios (lBrudereill2013h . 

DALI use s a r eacti on network de scribed in detail in 
IBruderer et al.l (l2012h and iBrudererl (l2013ll . It is based on a sub¬ 
set of the UMIST 2006 gas-phase network (I Woodall et al.l200^ . 
About no species and 1500 reactions are included. In addition 
to the gas-phase reactions, some basic grain-surface reactions 
(freeze-out, thermal and nonthermal sublimation and hydrogena- 
tion like g:0 —> g:OH —> g:H20 and H 2 /CH^ formation on PAHs 
(I.Ionkheid et al.ll2()()4 are included. The g:X notation refers to 
atoms and molecules on the grain surface. The photodissoci¬ 
ation rates are obtained from the wav elength-dependent cross- 
sections by Ivan Dishoeck et al.l (l2006l) . The adopted cosmic-ray 
ionization rate is {^2 - 5 x 10“*’ s“*. X-ray ionization (X-ray 
luminosity of 10^** erg s“*) and the effect of vibrationally exited 
H 2 are als o included in the netw ork. Elemental abundances are 
taken from lJonkheid et akl (l2006l) . 

3.3. Approach 
3.3.1. Physical model 

As a starting point for our models we ado pt the physical struc- 
ture suggested by [Andrews et al.l (l201 ill , as implemented by 
lBrudereir (l2013l) and app lied for a similar analysis of Oph IRS 
48 (IBruderer et all I 2 OI 4 . The gas surface density is assumed 
to be an exponential power-law following the time-dependent 
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Fig. 2. Generic surface density profile for the gas and dust. 

viscosity disk model v ~ W (iLvnden-Bell & Pringlel Il974t 
iHartmann et al.lll998h 

2(r) = Sf j exp 

and is defined by surface density 2^ at critical radius r^. The 
power-law index y is taken as 1, in line with the results for nor¬ 
mal disks with c ontinuous dust distri butions and other transi¬ 
tion disk studies (lAndrews et al.ll201 Ih. The previous g as study 
of observations in IRS 48 dBrudereret al.l 12014 ) demon¬ 
strated that the slope of the surface density profile could not be 
determined with their data, whereas the data in this new study 
have even lower S/N, so y is not varied. The outer radius is set 
to Tout. The dust density is taken to be S(r)/GDR, with GDR the 
gas-to-dust-ratio set to 100 initially. The inner dust disk is de¬ 
fined by the Z(r) profile, starting from the sublimation radius 
Gub = 0.07(L*/Lo)'^^ (assuming Tsub^lSOO K) out to rgap=l 
AU, scaled down by a factor (5dust to fit the near infrared part 
of the SED. The size of the inner disk is arbitrary as it does not 
have a significant effect on the near infrared excess. Between 
Tgap and Tcav, the cavity radius, the dust density is scaled down 
by a factor ^dustcav (the drop in dust density inside the cavity), 
while the gas density E(r) is scaled down by a factor (5gas. Figure 
12 shows the generic density structure and Table [3] explains each 
parameter. 

For the stellar photosphere a blackbody with the luminosity 
and temperature as in Table[T]is generated. Additional UV excess 
from accretion is added through a 7'acc=10 000 K blackbody for 
the Lace corresponding to the measured accretion rate (Table [1]), 
using 

Lacc(v) = (2) 

^^ acc 

with v) the Planck blackbody function, G the gravita¬ 

tional constant, M, and R, the mass and radius of the star, M the 
mass accretion rate and cr the Stefan-Boltzmann constant. For 


SR21, no accretion luminosity was added since only an upper 
limit for the accretion is measured. Note that observed values of 
M are uncertai n up to an order o f magnitude and often variable 
with time (e.g. lSalvk et alJl2013h . The interstellar radiation field 
incident on the disk surface is also included but does not affect 
the results presented here. 

The vertical structure follows a Gaussian distribution, de¬ 
fined by a scale height angle h — he . A disk in hydrostatic 
equilibrium will have ^-Q.25, while a flat {hjr - constant) disk 
will have ijj - Q. The near- and far-infrared excess in the SED 
relate directly with the resulting scale height h{r) at rsub and rcav, 
respectively. Our vertical structure is a simplification and does 
not account for hydrostatic equilibrium of gas or dust tempera¬ 
ture. One of the main effects of the hydro static equilibrium is a 
puffed up inner rim. Previous studies (e.g. [Andrews et aDl201lh 
often use a parametrized puffed up inner rim to fit the near- and 
far-infrared excess by adding a narrow scale height peak, but 
during our fitting procedure it became clear that a parametrized 
optically thick inner disk rim blocks a large amount of the UV 
irradiation, low ering the gas tempe rature significantly, as also 
demonstrated in IWoitke et al.l (l2009ll . 

For the dust op acities we use the same approach as 
[Andrews et al.l ([201 ih . consisting of a small (0.005-1 ;um) and 
large (0.005 - 1 mm) population of dust with a mass frac¬ 

tion /is in large grains. Settling is parametrized by defining the 
scale- height of the large grai ns as xK is. !i is reduced by a factor 
X (see [Andrews et all (1201 ll) for definitions). We initially fix fjs 
and X to 0.85 and 0.2 respectively, as adopted by [Andrews et al.l 
([201 ih . We include PAHs since many of our observed s ources 
show PAH features ([Brown et al.ll2007l [Geers et ^120061) but at 
a low abundance of 0.1% of the ISM abundance. The turbulent 
width is set to 0.2 km s '. Stellar parameters taken from the lit¬ 
erature are given in Table [1] 

3.3.2. Model fitting approach 

In the fitting procedure, we use an approach of manual fitting by 
eye, varying the parameters in a logical order, c onverging to a 
model that fits the data, similar to the approach in iBruderer et al.l 
(120141) . We do not use a x^ or Markov-Chain-Monte-Carlo 
(MCMC) method, as the computational time of the models is 
too long and the number of parameters too large. Although this 
approach reduces computational time, deriving uncertainties of 
model parameters, verifying the uniqueness of the fit, and esti¬ 
mating the correlation between parameters is not possible. The 
fitting parameters are defined in Table |3] 

1. Fit Sc and ^just roughly throu gh the SED and contin uum im¬ 
age, using the rcav value from lAndrews et aki (1201 ll) . 

2. Adjust Tsub, rgap and (Jdust if necessary for the NIR excess. 

3. Vary the scale height angle to fit the NIR and/or FIR excess 
in the SED, by changing he and if/. 

4. Vary rcav and re to match the null and slope of the submil¬ 
limeter dust visibility curve. 

5. If the MIR/FIR excess is too high, adjust the size distribution 
by increasing /s and x- 

6. Find the highest possible value for i5dustcav (drop of dust den¬ 
sity inside the cavity) that is still consistent with the dust 
continuum. 

7. Check the fit to the CO intensity map, spectrum and visibil¬ 
ities assuming the gas-to-dust ratio (GDR) to be 100 in the 
outer disk. If necessary, truncate the outer radius to fit the 
spectrum. 
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Table 3. Parameters fitting procedure 


Category 

Parameter 

Value 

Description 

Surface density profile (E(r)) 

Fc 

free 

Critical radius 



free 

Gas density at r^ 


y 

1 

Surface density gradient 


GDR 

100 

Gas to dust ratio 

Radial structure 

^sub 

fixed by formula 

Sublimation radius 


rgap 

free 

Size inner disk 


^cav 

free 

Cavity radius 


out 

free 

Outer radius 

Scaling E(r) 

^dust 

free 

Dust density drop inner disk (rjub-rgap) 


^dustcav 

free 

Dust density drop inside cavity (rgap-rcav) 


^sas 

free 

Gas density drop inside cavity (rsub-rcav) 

Vertical structure (h{r)) 

he 

free 

Scale height angle at r^ 



free 

Flaring angle 

Dust properties 

sg-pop 

0.005-1 yum 

Small dust grain population 


Ig-pop 

0.005 pm-l mm 

Large dust grain population 


h 

free 

Fraction large grains 


X 

0.2 

Scale height large grain fraction 


PAH 

0.1% 

PAH fraction w.r.t. ISM abundance 

Ionization 

Cm 

5xl0-''s-‘ 

Cosmic ray ionization rate 


X-ray 

10^** erg s“' 

X-ray ionization rate 


8. Explore small variations in the vertical parameters to change 
the temperature structure, required to fit the optically thick 
'^CO emission. 

9. Finally, keep all other parameters constant while varying dgas 
to find the lowest possible value for 5gas that is still consistent 
with the CO data. 

We stress that this approach does not result in an unique 
model, due to the degeneracies in parameter choices and due 
to '^CO being optically thick. However, we consider this ap¬ 
proach to be reasonable considering the different data sets taken 
into ac count, where ALMA is the crucial one. ICarmona et al.l 
(l2014l) also demonstrated a modeling procedure combining dif¬ 
ferent datasets. The submillimeter continuum sets firm con¬ 
straints on the total flux, the cavity size (the null in the visibil¬ 
ity )_and_the_di^exteiU (through Fc) of the disk, as also shown 
in [Andrews et al.l (1201 ih . The mid infrared part of the SED and 
the are both sensitive to the combination of temperature 
and density structure of the disk, rather than just the density, as 
both of these emission regimes are mostly optically thick. The 
temperature is controlled by the adopted vertical structure in the 
disk, as this determines the direct irradiation by the star, espe¬ 
cially at the cavity wall. Therefore, several iterations were re¬ 
quired to find a compromise for a decent fit of both the SED and 
i^CO. 

3.3.3. CO emission inside the cavity 

As the CO emission inside the cavity only partly depends on den¬ 
sity, the derived bgas parameter is merely a constraint within the 
derived physical model. As discussed above, there are several ef¬ 
fects that influence the strength of the emission: most important 
are gas temperature and CO abundance which are very sensitive 
to the small dust grains inside the ca vity. Different ef fects can 
play a role, as discussed extensively in lBruder^ (l2013h . 

- Although the bulk of the '^CO emission is optically thick, the 
drop in gas surface density can still be constrained to within 
an order of magnitude with spatially resolved observations. 
Optically thick emission will slightly drop with decreasing 


bgas because the r = 1 surface moves deeper into the disk, 
so the emission traces colder regions in the disk, as d emon- 
strated in Section 5.2 and Fig. 9-10 in iBruder^ (|20 1 3h . 

- An inner disk will cast a shadow on the gas inside the 
cavity, depending on the amount of dust and scale height. 
This will lower the gas temperature inside the cavity. On 
the other hand, the inner disk shields the stellar UV radia¬ 
tion that would photodissociate CO, so the CO abundance 
is increased. The structure of the inner disk is not well con¬ 
strained from the available data, while it may have significant 
effects on the heating of the outer disk. 

- The presence of small dust grains and PAHs inside the cav¬ 
ity shields stellar UV radiation and decreases the photodis¬ 
sociation, thus increasing the CO abundance. The dust con¬ 
tributes to the heating and cooling of the gas, depending on 
the density. For high gas density, the gas temperature de¬ 
creases when dust is added due to cooling by gas-grain col¬ 
lisions. For low gas density, the gas temperature decouples 
and becomes higher than the dust temperature due to photo¬ 
electric heating on grains and PAHs. 

- The increase of the gas temperature (increase of small dust 
grains) or an increase of the UV radiation (decrease of small 
dust grains) will also lead to an increase of vibrationally 
pumped H 2 , again increasing the formation of CO. 

- The presence of small dust grains and PAHs increases the 
formation of CH^ on the grain surfaces, which leads to in¬ 
creased formation of CO. 

The consequences for the CO emission of adding small dust, ei¬ 
ther inside the entire cavity or just in the inner disk, strongly de¬ 
pends on the details of the structure of the disk. This is the reason 
why also the SED and submillimeter visibilities are included in 
the modeling procedure, and why a full physical-chemical model 
is required for analysis of the CO emission. 

4. Results 

Our best fitting results for the dust and gas density structure are 
presented in Figures[3]to|5]and Table|4] In Figure[3]the output of 
the physical model is compared to the SED, the 690 GHz con- 
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Fig. 3. Modeling results and observations of the dust surface density, comparing bdustcav ranging between 10^^ and as indicated 
in the right panel. The observations are plotted in black. Left: Spectral Energy Distribution; middle left: Amplitude of the 690 GHz 
continuum visibility for the deprojeeted baselines. The null line is indicated with a dashed line; middle right: Normalized intensity 
cuts through the major (bottom) and minor (top) axis of the 690 GHz continuum image. The model images are convolved with the 
same beam as the ALMA observations; right: The dust surface density profile. Indicated are the ddust, the drop in density to fit the 
inner disk through the near infrared emission, and bdustcav, the minimum drop in dust density inside the cavity to fit the observations. 
In SR21 (top panel) the region inside 7 AU is assumed to be empty due to additional constraints (see main text). In all other disks 
(5dustcav is found to be at most 10^'*. All other disk parameters are as listed in Table|4] 
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Fig. 4. Modeling results and observations of the '^CO emission for the best fitting physical model derived from the dust (see Figure 
[3 and Table m, with bgas=l (no drop in gas surface density inside the cavity). The observations are black, the models in purple. 
Left: '^CO spectrum integrated over the entire disk; Center: Intensity cuts through the major (bottom) and minor (top) axis of the 
zero moment map. The noise level is indicated by the gray zone. The model images are convolved with the same beam as the 
ALMA observations. ; Right: Amplitude of the integrated *^CO visibility for the deprojected baselines. The null line is indicated 
with a dashed line and error bars indicate the noise. 


tinuum visibility curve and the normalized major/minor axis cuts 
of the 690 GHz continuum image (left panel of Figure [T]i. Data 
poi nts of the SEP have been dere ddened using the extinction law 
by IWeingartner & Draind (1200 ih with Rv=5.5. The SED con¬ 
strains the vertical structure of the disk, whereas the visibility 
curve determines the cavity size and extent through the location 
of the null. The image cuts show the direct comparison of the 
image with model. Within each figure, the bdustcav value (drop 
of dust density inside the cavity) is varied with values between 
10“^ and 10 to constrain the minimum drop in dust density 


inside the cavity. In the right most panel the density structure of 
the model is presented. 

Reasonable fits are found for the continuum visibilities and 
SED within the constraints of our physical model for all targets. 
The inferred disk masses a re comparable to prev ious findings 
within a factor of a few, e.g. [Andrews et al.l (1201 Ih . The models 
with bdustcav=10“^ overproduce the emission inside the cavity in 
all cases, as seen both in the mid infrared part of the SED and 
in the dust continuum image. Therefore, ^dustcav is constrained to 
be at most a factor of 10“"^, and even as low as 10“® for J1604- 
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Table 4. Results fitting procedure for the gas density profile of each transition disk. 


Target 

fc 

(AU) 

2. 

(g cm-2) 

^dust 

(lO-'Mo) 

-^gas 

(lO-'Mc,) 

he 

(rad) 


<^dust 

f cav 

(AU) 

^ sub 

(AU) 

rgap 

(AU) 

/i. 

rout 

(AU) 

^gas 

*5dustcav 

SR21 

15 

400 

0.12 

12 

0.07 

0.15 

10-'’ 

25 

0.18 

1 

0.85 

100 

10-2 

10-5 

HD135344B 

25 

200 

0.17 

24 

0.15 

0.05 

10-2 

40 

0.18 

0.25 

0.95 

150 

10-' 

< lo-'* 

LkCal5 

85 

34 

0.95 

103 

0.06 

0.04 

10-5 

45 

0.08 

1 

0.98 

400 

10' 

< 10-* 

RXJ1615-3255 

115 

60 

3.5 

470 

0.04 

0.2 

10-'’ 

20 

0.08 

1 

0.85 

200 

> lO-** 

< 10-5 

SR24S 

15 

1200 

0.35 

- 

0.12 

0.01 

10-^ 

25 

0.14 

1 

0.98 

- 

- 

< 10-* 

J1604-2130 

60 

12 

0.095 

20 

0.065 

0.68 

10-‘ 

70 

0.04 

0.06 

0.98 

400 

10-5“ 

< 10'’ 


Notes. This value refers to the drop inside the inner cavity of 30 AU (dgas 2 in the text). 


2130. In the modeling of the gas, we set (5dustcav to zero, except 
for SR21 where ddustcav is constrained by the visibilities. For the 
gas modeling of HD135344B and J1604-2130, we include some 
additional dust between and 30 AU as shown in the dust 
density profile in Figure |5] 

For the analysis of the gas density, we compare the result 
of the adopted disk model to the '^CO spectrum, the intensity 
cuts through the zero moment map and the visibility curve 
in Figure H] SR24S is not included in this part of the analysis 
because the CO emission is too heavily confused by large scale 
emission (Fig. lA.l| i. The models are calculated assuming (igas=l, 
i.e. no drop in the gas surface density inside the dust cavity, ex¬ 
cept for J1604-2130, which shows a resolved drop in the CO 
emission inside 30 AU, which is adopted already. Because the 
CO visibilities do not show a clear null and the disks are not 
fully symmetric in the gas, the visibility curve was not used to 
fit our data exactly; it is merely included here to show that the 
emission from the moment map is well retrieved from the visi¬ 
bilities. 

The models reproduce the CO emission well in the outer part 
of the disk (r > r^^y). There was no need to change the gas to 
dust ratio of 100 in any of the cases in the outer disk. This con¬ 
clusion needs to be confirmed by optically thin CO isotopologue 
observations. The images and visibilities show that the outer ex¬ 
tent of the dust and CO emission can indeed be fit with the same 
physical model, although better S /N observations are required 
to confirm this. 

For SR21, HD135344B and LkCalS, the initial model over¬ 
produces the central part of the intensity cut by a factor 1.5-2, 
suggesting a decrease of gas surface density inside the dust cav¬ 
ity. In Figure |5] we present the same physical model with dgas 
varying from 10° to 10^“* in steps of 10 to compare with the ob¬ 
servations. The emission drops by less than a factor of 2 in each 
of these steps due to the optical depth of the '^CO line, but within 
the noise level of the moment map this is sufficient to constrain 
<5gas to within an order of magnitude due to the different gas tem¬ 
peratures at the T = 1 surface in the disk. This result is different 
from Oph IRS 48, where the *^CO line wings were found to be 
optically thin (iBruderer et alJ20l4l) . Only the emission in J1604- 
2130 inside 30 AU radius becomes optically thin. The range of 
i5gas is indicated in the right panel of Figure |5] and in Table |4] 

In the following section we discuss the properties of each 
disk in more detail. 

4.1. Results of individual targets 
SR21 

It is possible to fit the SED and amplitudes of the continuum 
visibilities simultaneously at least up to baselines of 600 kd. The 
small inconsistency at longer baselines is possibly related to the 


dust asymmetry (iPerez et al.ll2014l) which is not taken into ac¬ 
count. In order to fit the SED, we use for this disk an alternative 
dust opacity table for the small grains, without small silicates, 
because the normal dust table results in a strong silicate feature 
that is not seen in the data. We find a cavity size of at most 25 
AU, which is significantly smaller than the v alue of 36 AU pre¬ 
viously found in the SMA data at 345 GHz (iBrown et alJl200^ 
[Andrews et al.ll2(Hl]) . It remains unclear whether this is due to 
differences in the modeling approach or a physical effect, since 
the 345 GHz emiss ion traces somewhat l arger dust grains (see 
also 5.2). Similar to [Andrews et akl ([201 Ih we find that this disk 
has a remarkably small radius: is only 15 AU, resulting in no 
detectable emission outside of 75 AU radius. 

Both hot CO line emission ([Pontonnidan et al] [2008h and 
mid infrared interferometry (Benisty, private comm.) suggest a 
hot ring at 7 AU and no material inside. Also s cattered light 
show s that there is small dust down to <14 AU ([pollette et al.[ 
[2013[) . Therefore, we assume that there is no gas between the star 
and 7 AU for SR21 (but still an inner dust disk due to the NIR 
excess), and dust with only a small grain population between 7 
AU and rcav The dust between 7 AU and the dust hole cavity of 
25 AU is partly constrained by the visibility curve: without the 
additional dust the second null is much further out, although this 
may also be due to the asymmetry. In constraining ^dustcav, this 
drop is assumed to start from 7 AU rather than from rgap (taken 
as I AU). 

The central '^CO emission in the intensity cut is overpro¬ 
duced by about a factor of 2 in the initial model, especially inside 
the cavity region (Pig. |4]i. Only a deep drop in 5gas of 2 orders 
of magnitude can match the peak of the zero moment map, al¬ 
though the intensity cut is somewhat wider than the data. A mi¬ 
nor dip of the 7 AU gap is seen in the intensity cut of the model 
integrated '^CO emission; the gap is marginally resolved, but 
the dip is within the error bars of the observations so it cannot be 
confirmed with these data. 

HD135344B 

The continuum visibility is reasonably well fitted with a model 
with a 40 AU cavity, slightly smaller than the 45 AU cavity found 
in previous wo rk with 345 GHz SMA data ([Brown et al.[[200^ 
[Andrews et al.[l2()TTI: [Carmona et al.[ 120141) . An increase of the 
large grain fraction and large grain scale height fx - 0.8) is 
needed to fit the peak of the far infrared part of the SED. 

We find a smaller critical radius of only 25 AU and thus 
smaller extent of the disk compared to [Andrews et al.[ ([201 ih 
(who found rc-55 AU), but the ALMA data quality is much bet¬ 
ter than the SMA data. Another main difference in the modeling 
procedure is the different vertical structure of this disk in our 
parametrization (no puffed up inner rims), although this is not 
expected to change the emission in the outer disk. Scattered light 
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Fig. 5. Modeling results and observations of the *^CO emission for the best fitting physical model for gas and dust (Figure HI 
comparing (5gas (gas density drop inside the dust cavity), dgas is varied as 10“, ... , lO^'* as blue, light-blue, green, yellow and 

red in the left and middle plots, while the observational data is in black. Left: The spectrum integrated over the disk area. 
Middle: The intensity cuts through the major (bottom) and minor (top) axis of the '^CO zero moment map for different bgas. The 
noise level is indicated by the gray area. The model images are convolved with the same beam as the ALMA observations. Right: 
The best fitting result for possible values of 5gas and ddustcav in red and blue respectively. 


images show that there is small dust inside the millimeter dust 
cavity do wn to 30 AU for HD 135344B (millimeter dust cavity is 
~46 AU) (iGarufi et al.ll2013l) . Therefore we assume a finite drop 
in dust density between rcav and 30 AU. The scattered light im¬ 
ages, SED and visibilities do not constr a in the amo unt of dust: 
we choose a drop of 10~^. iGarufi et alJ (1201 3h and iMuto et alJ 
(1201 2h show signatures of spiral arms in the scattered light im¬ 


ages. The ALMA data have insufficient resolution to detect these 
arms if they were visible at millimeter wavelengths. 

The initial model for this disk overproduces the CO emission 
inside the cavity. A dgas of 10“' - 10“^ is sufficient to fit the 
peak in the moment map. A decrease in gas density inside the 
cavity was previously inferred from SM A CO line obseryations 
(iLyo et alJ201 Ih . lCarmona et al.l (12014l) require a surface density 
that increases with radius inside the cayity to fit the royibrational 
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CO lines in this disk. Fitting these lines is beyond the scope of 
this paper, but the presence of a drop in surface density may be 
similar to this increasing density profile, as discussed in Section 
5.1. Motivated by the dust structure we also try a gas model with 
a cavity size of 30 AU rather than 40 AU. It turns out that there 
is no significant difference in the outcome for 30 AU or 40 AU 
gas cavity radius. As Figure|6]shows, (5gas = 10 * - 10 ^ also fits 
best for a 30 AU gap. 



Cut (") 


Fig. 6. Modeling results and observations of the emission 
for the best fitting physical model for HD135344B for a 30 AU 
rather than 40 AU cavity for gas and dust (Figure|4]i for different 
values of (5gas (gas density drop inside the dust cavity), dgas is 
varied as 10*’, 10“', ... , lO^'' as blue, light-blue, green, yellow 
and red. The observations are plotted in black. The panel shows 
the intensity cuts through the major (bottom) and minor (top) 
axis of the '^CO zero moment map. 


dust, as seen in recent ALMA data of HD 100546 (IWalsh et al.l 
12014 . 



Cut (”) 


Fig. 7. Modeling results and observations of the 690 GHz con¬ 
tinuum of RXJ1615-3255. The panel shows the intensity cuts 
through the major (bottom) and minor (top) axis of the contin¬ 
uum image, zoomed in on the lowest flux levels (< IQcr). The 
observations are plotted in black with the Icr errors in gray. In 
red and blue the models with and without a dust gap between 110 
and 130 AU are presented. The dotted and dashed lines indicates 
the zero flux and 3cr limits, respectively. 


The model with 5gas=l does not reproduce the high velocity 
emission seen in the CO spectrum, but the S/N is low so the 
significance of this remains debatable. The intensity cuts show 
that the gas density inside the dust hole can be up to 4 orders of 
magnitude lower before the CO becomes optically thin. 


LkCalS 

LkCal 5 is a well studied transition disk with a large cavity of 
50 AU (llsella et al.ll2012h . We model LkCal5 with a rather flat 
but massive disk structure. Our derived disk dust mass is 3 times 
higher than typically found in the literature, and the 690 GHz 
continuum is indeed optically thick, which explains the apparent 
asymmetry in the continuum image: the north west side of the 
disk is brighter than the south east due to the geometry. 

Just like HD135344B, the initial model for this disk overpro¬ 
duces the CO emission inside the cavity, indicating a dgas of at 
least 10“'. 

RXJ1615-3255 

RXJ1615-3255 does not reveal a cavity in the continuum image, 
and also the visibilities show only a hint of a null at the longer 
baselines (>500 kd). Our best fitting model has a 20 AU cav ity, 
somewhat smaller than the best fit in [Andrews et al.l (1201 ll) of 
30 AU. A model without cavity can not fit the SED which shows 
a clear indication of a deficit of dust close to the star. This in¬ 
dicates that the Band 9 continuum might be optically thick. An 
interesting aspect of this target is the outer part of the disk: a 
narrow radial dip is seen in the image at 0.6” from the center 
(corresponding to 110 AU at the distance of this disk). Figure 
|7] shows a zoomed-in image of the intensity cut from Figure |3] 
This hints at the presence of a dust gap in the outer part of the 
disk between 110 and 130 AU, or an additional outer ring of 


SR24S 

SR24S is an inclined disk with a just resolved cavity of about 30 
AU. The long wavelength part of the Spitzer IRS spectrum and 
the MlPS-24 flux are k nown to be confused b y its nearby com¬ 
panion SR24N (see e.g. [Andrews et al.1l201 ll) . thus these points 
are only used as upper limits in our fit. The visibility curve be¬ 
comes inconsistent with our model for baselines >350 kd, hint¬ 
ing at substructure inside the cavity. As we do not have any ad¬ 
ditional information on dust inside the cavity of this disk we do 
not explore this further. 

J1604 

J1604-2130 is a well-studied trans ition disk at 34 5 GHz, both 
with SMA ([Mathews et al.l 120121) and ALMA (IZhang et al.l 
I20l4 . Both authors prov i de ph ysical models for both the dust 
and the gas. IZhang et"^ (1201 4 uses the same ALMA dataset 
discussed here, but does not use a full physical-chemical model 
for the gas analysis. Also, in both studies the near infrared excess 
is not taken into account. The issue is the apparent variability: 
whereas the Spitzer IRAC photometry (taken in 2006 and used 
by Mathews) shows no near infrared excess, the Spitzer IRS and 
WISE photometry (taken in 2007 and 2012, respectively) do re¬ 
veal the excess. Zhang et al. note this discrepancy and state that 
the near infrared excess can be fit with an optically thick inner 
disk ring, but do not include this in their model. Since the inner 
disk has a strong influence on the CO emission inside the cav¬ 
ity (and on the shadowing of the dust cavity wall), we choose to 
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include the more recent photometry with near infrared excess in 
our SEP and fit the model s accordingly. 

iMathews et al.l (1201 2l) introduce a two drop dust density 
model to fit their SMA data which is also adopted bv IZhang et alJ 
(1201 4l) in the analysis of the ALMA data: one drop at 75 and one 
at 20 AU. We find that such a two-step decrease is indeed re¬ 
quired to fit the continuum visibilities and adopt the same struc¬ 
ture, although we set the inner drop radius at 30 AU to be con¬ 
sistent with the gas cavity. On the other hand, we require a much 
larger critical radius of 60 AU rather than 10-20 AU as used by 
Mathews and Zhang. This is likely the result of shadowing due 
to the presence of an inner disk. Our disk dust mass of 0.1 Mjup 
is comparable with their estimates. 

In the gas a clear signature for a drop at 30 AU is see n in the 
CO moment map, which was modeled by IZhang et~an (l2014l) 
in a parametrized gas model of a power-law temperature profile 
with an empty gap. Their result is consistent with that of our full 
physical-chemical model and the derived surface density profile 
is similar (right panel of Figure©. We can constrain the inner 
density drop to (5gas2 = 10“^ (see Figure©, although this depends 
on the inner cavity size: a smaller cavity size would require an 
even deeper drop, a larger size a less deep drop. On the other 
hand, 5gas (the drop inside the dust cavity radius) is 1, there is no 
visible drop in density between 30 and 70 AU. 




Cut (") Radius (AU) 

Fig. 8. Left: Modeling results and observations of the *^CO 
emission for the best fitting physical gas model for J1604-2130 
for the inner 30 AU cavity comparing (5gas2 (gas density drop 
inside the 30 AU dust cavity). (5gas2 is varied as 10*^, 10“', ... , 
10“® as blue, light-blue, cyan, green, yellow,orange and red (bgas 
is taken as 1). The observations are plotted in black. The panel 
shows the intensity cuts through the major (bottom) and minor 
(top) axis of the *^CO zero moment map. Right: Comparison 
betwee n the density profil e derived in this study and the profile 
used in IZhang et al.l (120141) . 


5. Discussion 

5.1. Implications of gas density drop 

The key parameter in our analysis is bgas, the drop in gas density 
inside the dust hole, in comparison with bdustcav For each target, 
the minimum drop in dust density bdustcav is significantly lower 
than the possible values of bgas. Figure ©and Table ©show the 
range of possible values for these parameters for each target. 

The significantly larger drop in dust density compared with 
the gas is consistent with the dust trapping scenario by a com- 
panion, with gas p resent inside the millimeter-sized dust cavity 
(iPinilla et al.ll2012h . For the sources studied here, photoevapora¬ 
tion by itself is unlikely based on just the large cavity sizes and 


high accretion rates (lOwen et al.ll201 it iRosotti et al.ll2013l) and 
our results further strengthen this hypothesis. Especially SR21, 
HD135344B and LkCal5 require a minimum drop in the gas 
surface density of a factor 10-100 within the gap and are thus 
promising candidates for disks with embedded planets. The drop 
is deep enough to rule out the dead zone scenario, which has a 
drop that is much less than an order of magnitude (iLvra et al.l 
l2015h . On the other hand, one has to be careful with overinter¬ 
pretation at this point: the CO emission inside the gap is sensitive 
to the amount of small dust inside the cavity due to its shielding 
effect, which can also explain the apparent drop in CO emis¬ 
sion. We have explored this within the constraints of our physi¬ 
cal model, but any addition of dust or PAHs inside the cavity that 
is sufficient to lower the CO emission at the same time increases 
the mid infared continuum emission in the SED to levels that 
are inconsistent with the data. The combination of SED mod¬ 
eling and spatially resolved *^CO emission thus provide unique 
constraints on the temperature and small dust grain composition. 
On the other hand, the vertical structure also influences the tem¬ 
perature structure and thus the *^CO emission and mid infrared 
continuum emission, and this degeneracy remains uncertain. 

The depth of the gas density drops at the cavity radius is 
still modest, only a factor bgas ~ 10 * - 10“^. Comparison 
with the plane t-disk interaction m odels of Jupiter mass planets 
in Eigure 1 in iPinilla et al.l (12012h indicates that any embedded 
planets responsible for this drop are unlikely to be more massive 
than 1-2 Jupiter masses, and even then only with high viscosity 
(a ~ 10“^). A high visco sity is inconsistent with the pres e nce of 
long-lived vortices (e.g. lAtaiee et al.ll2013l) . iPerez et ^ (1201 4l) 
interpret the asymmetries of SR21 and HD135344B as possi¬ 
ble vortices, by comparing their analytic solution f or the asym¬ 
metric dust structures with the vortex prescription of iLvra & LinI 
(l2013h . The lack of a dee p gap in the CO argues against this in¬ 
terpretation. [pinniaer^jloT^ also conclude that for SR 21 the 
vortex scenario is unlikely by analyzing the observed dust distri¬ 
bution at different wavelengths with hydrodynamic and dust evo¬ 
lution models. However, a smaller gas cavity size with a deeper 
density drop would also fit our observations and be consistent 
with low viscosities. Eurther high spatial resolution observations 
of CO isotopologues are needed to disentangle cavity size and 
gas density drop. 

The planet-disk interaction models actually predict rounded 
off gradients rather than steep vertical drops such as 
parametrized in our model. A rounded-off dust wall was indeed 
inferred from mid in frared interferometry data of HD 100546 
(iMulders et al.ll201^ by fitting the visibility curve, but the S/N 
of our data is insufficient to test this. One important conclu¬ 
sion of IPinilla et al.l (1201 2h is that the gas cavity radius (loca¬ 
tion of the planet) is expected to be up to a factor of 2 smaller 
than the dust cavity radius (location of the radial pressure bump 
where the dust is trapped ). This was indeed obse rved for J1604- 
2130 (IZhang et al.l2014i this studvl. Oph IRS 48 (iBruderer et al.l 
l2014t) and HD 142527 (iPerez et al.|[2015h . As stated above, our 
estimates of the drop in gas density inside the cavity would likely 
increase if we would assume a smaller gas cavity size, but the 
spatial resolution of our data is insufficient to constrain this. 
Eor HD135344B we show the outcome for a gas cavity of 30 
rather than 40 AU (Pig. ©, but for this radius the bgas is still 
limited to 10 * - 10“^. Eor J1604-2130, we find clear evidence 
for an inner drop at 30 AU, and no additional dr op between 30 
and 7 0 AU of 10“*, as found for Oph IRS 48 (IBruderer et al.l 
l2014l) . A structure with a double drop or smaller gas cavity size 
hints at the presence of more than one companio n, or alterna- 
tively a shallow increasing slope such as seen in IPinilla et al.l 
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Table 5. Location of the null in the visi bilities of models an d 
SMA observations at various frequencies (lAndrews et al.ll2()Tlll . 


Target 

Model-690 
null (U) 

Model-345 
null (U) 

Observed-345 
null (kd) 

SR21 

252 

247 

200-220 

HD135344B 

196 

204 

160-180 

LkCal5 

171 

155 

140-160 

RXJ1615-3255 

484 

349 

250-300 

SR24S 

255 

260 

240-260 


(1201 2l l2015h . This possibility can not be excluded for any of 
the other sources in t his study with the available observations. 
ICarmona et al.l (l2014t) suggest a gas surface profile increasing 
with radius inside the cavity to match the gas and dust density 
distribution of HD 135344B, although this is only constrained by 
observations for the inner few AU. Modeling the shape of the 
gap is beyond the scope of this paper. Spatially resolved obser¬ 
vations of CO isotopologues will provide better constraints on 
the density profile. 

5.2. Cavity size and dust distribution 

Another interesting aspect of our results is that all disks have sig¬ 
nificantly smaller cavity sizes (by 5 -10 AU) in the 690 GH z con¬ 
tinuum data than those derived by [Andrews et al.l (1201 ih based 
on the 345 GHz continuum observations taken with the SMA. 
The cavity size as derived with the ALMA continuum visibili¬ 
ties (location of the null) is constrained to better than +5 AU. 
The S /N of the SMA data is much lower, so the cavity size is 
less precisely constrained but the systematic trend is significant. 
For comparison, we list the location of the first null for both the 
690 GHz and 345 GHz continuum of our adopted model in Table 
|5]and compare those with the null in [Andrews et al.l (1201 ll) . 

Note that the null as derived from the models can shift 
slightly for the two frequencies, due to difference in dust opacity 
at the observed wavelengths. The null in the observed 345 GHz 
continuum is significantly smaller for SR21, HD135344B and 
RXJ1615-3255 than that at 690 GHz, indicating a larger cavity 
size at this longer wavelength. This hints at a radial variation of 
the dust size distribution, as the continuum emission is sensitive 
to larger grains at longer wavelengths. A more concentrated dust 
ring (i.e. a narrower ring and thus a larger dust cavity size) is con¬ 
sistent with the dust trapping scenario, because larger dust grains 
are further decoupled fr om the gas and thu s more affected by ra¬ 
dial drift and trapping (iBrauer et al.ll200^ . Higher S/N ALMA 
Band 7 observations are required to confirm this. 

6. Conclusions 

In this work, we have analyzed high spatial resolution ALMA 
submillimeter observations of '^CO line emission from 6 transi¬ 
tional disks using a physical-chemical model. By comparing the 
SED, the 690 GHz continuum visibilities and the *^CO emission 
simultaneously, we derive a physical model of the gas and dust. 
With the model we can set constraints on the surface density pro¬ 
file of the dust and gas and specifically, the amount of dust and 
gas in the cavity. 

1. All disks show clear evidence for gas inside the dust cavity. 

2. The gas and dust observations can be fit with a surface den¬ 
sity profile with a steep density drop at the cavity radius, 
taking a gas-to-dust ratio of 100 in the outer disk. 


3. The combination of SED and spatially resolved '^CO inten¬ 
sity fitting sets constraints on the vertical and thus tempera¬ 
ture structure of each disk. 

4. All disks except SR24S have a potential drop of 1 or 2 orders 
of magnitude in the gas surface density inside the mm-sized 
dust cavity, while the minimal drop in dust surface density 
inside the cavity is at least 3 orders of magnitude. 

5. J1604-2130 has a deep resolved gas cavity that is smaller 
than the dust cavity. Eor the other disks it is possible that the 
gas cavity radius is smaller than the dust cavity radius, in 
which case the density drop will be deeper, but this can not 
be constrained with the available data. 

6. The derived density profiles suggest the clearing of the cavity 
by one or more companions, trapping the millimeter dust at 
the edge of the cavity. 

7. Our model for J1604-2130 i s mostly con s istent with the pro¬ 
posed physical structure by [Zhang et al.l (120141) . derived us¬ 
ing a parametrized temperature model. 

8. The continuum of RXJ1615-3255 shows an additional dust 
ring around 120 AU. 

9. The derived cavity sizes of the millimeter dust (and loca¬ 
tion of the null in the visibilities) are consistently smaller for 
the 690 GHz/0.44 mm continuum than for the 345 GHz/0.88 
mm continuum. This is consistent with the dust trapping sce¬ 
nario, because trapping is more efficient for larger dust grains 
probed at longer wavelengths. 

The derived physical model can be used as a start for the 
analysis of CO isotopologue observations, putting better con¬ 
straints on the gas density inside and outside the cavity. We have 
recently obtained CO isotopologue ALMA observations for two 
disks of the sample in this study (van der Marel et al. in prep.). 
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Table A.l. Velocity range within the cavity radius 


Target 

Velocity (km s“‘) 

SR21 

<1.25 and >4.0 

HD135344B 

<5.0 and >9.0 

LkCalS 

<2.2 and >10 

SR24S 

<1.25 and >8.25 

RXJ1615-3255 

<-0.3 and >9.5 

J1604-2130 

<4.1 and >5.3 


SR21 



1.0 0.5 0.0 -0.5-1.0 
ARA (") 


HD135344B 



1.0 0.5 0.0 -0.5-1.0 
ARA (”) 


V (km s ') 

- 1.0 1.0 2.9 4.9 6.9 




1.0 0.5 0.0 - 0.5 - 1.0 

ARA (") 


LkCalS 



2 1 0 - 1-2 
ARA (") 


V (km s 

2.0 4.1 6.2 B .3 10.4 


1.5 
1.0 
^ 0.5 

O 0.0 
0 ) 

Q 

< - 0.5 



1.5 1.0 0.5 0 . 0 - 0 . 5 - 1 . 0 - 1.5 
ARA (”) 


Fig. A.l. *^CO channel maps for each observed target. Overlaid in white contours are the Keplerian velocity prohles for the derived 
inclination and given stellar mass. The yellow ellipse indicates the dust hole radius and the yellow contours the 5cr outer radius of 
the dust continuum. 
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Fig. A.l. *^CO channel maps for each observed target (continued) 
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